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Abstract 

Using numerical modeling investigated interaction of solitary waves (solitons) of 
the regularized long wave equation. For reception the stable model of the nonlinear 
medium are used methods of the linear prediction and progressive approximation. 
By modeling was determined that depending on ratio of velocities of the solitons and 
the form of highest derivatives balance is possible self-organization of the medium 
nonequilibrium state as formation of shock waves and stable on the form solitary 
waves, created as a result of full or partial mutual penetration of the solitons. Is 
possible also aggregation of the solitons in third wave. The shock waves can pass 
into other possible resonance state as a wave front with stable amplitude, which 
precedes developing in singularity negative front. 

Spatially separated solitary waves (SW) are in balance with nonlinear medium if they 
are resonance to it. At rapprochement of the SW the balance is saved when multiwave 
field is resonance too. Otherwise the wave field dynamic is difficult predictable. An in- 
teraction of the solitary waves of the nonresonance wave field is possible to analyse by 
numerical modeling. As analyse object we shall choose nonlinear medium, described by 
the regularized long wave (RLW) equation 

u t + iu x - au xxt + 2f3uu x = 0, (1) 

where u = u(x,t), 7 = 1, subscript symbols mark the partial derivatives on spatial x and 
temporal t variables. The RLW equation describes haracteristic of the nonlinear medium 
in more broad dynamic range in contrast with Korteveg - de Vries's (KdV) equation, 
however, it is't solvable by the inverse scattering method and other algebraic methods and 
so has't A-solitons solutions. The RLW equation is considered in literature with pairs of 
coefficients: (a; (3) = (1; 0.5) pQ, (l;-3) j2] and (0.5;0.125) [Sj. Substitution u = w x and 
Hopf - Cole's transformation t lJ as w = X(\nT) x , where T(x,t) = exp(a(x—vt)) + l, A, a — 
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constants, v - velocity, yield the condition of the highest derivatives balance GaT t +\(3T x = 
0, from which follows that the velocity 

v = \p/6a. (2) 



From linear part JIJ follows that spectral parameter a = \/ (v — j)Jav. The solution of 
the RLW equation with account of the expression for velocity is 



3 2 / 1 l V -7, 



u(x,t) = —(v -7)sech l-^—-( x -vt)j. (3) 

When v > 7 the function u(x,t) describes the solitary wave, which on nature of the 
dependencies of the amplitude from velocity, locality and conservation of the form at 
spatio-temporal translations corresponds to the determination of the soliton [Tj. When 
v < 7 the solution 

is periodically singular. We shall research an interaction of the one-soliton solutions @ 
for different ratios of the velocities and pairs of the coefficients (a, (3) using numerical 
simulation of the dynamic of an initial layer, defined as two waves, given on a square grid 
with intervals of sampling Ax, At on space and time accordingly. 

The numerical modeling of nonlinear processes is characterized ambiguous results, 
since presence power-mode function implies minimum two decisions. If chosen decision will 
not matched on the dynamic with previous, that the numerical scheme can lose stability. 
Since a nonlinearity reveals itself as influence of a sought value on itself, that appears 
the problem of the choice of the supporting area for generation of following values of the 
net under influence not yet defined value. Therefore computing model of the nonlinear 
process must be built on scheme of the progressive approximation with restriction on area 
of possible result values. 

Invariance of the solitons to spatio-temporal translation is expressed by equation 

u t + vu x = (5) 

with the dispersion equation 00 — v k = 0. Let Uij = u(jAx,iAt) denotes the samples of 
the discrete wave field, where i , j = 0,1,..., and Vq = Ax/ At is the velocity of the moving 
wave, matched with the samples net. As follows from equation (JSJ), for the soliton with 
velocity v = mt> , where m - integer number (the velocity index), is obviously identity 
By applying ^-images of the shift on time z t and space z x instead of finite 
differences the discrete analogue of the equation © can be written as 1 — Ztz x m = 0. For 
two solitons, moving with multiple t>o velocity v\ = ra\VQ and V2 = rri2Vo, the z-image of 
the dispersion equation {oj— Vik)(u— v 2 k) = transforms to 1— ZfZ™ 1 — ZtZ™ 2 + z\ ^™ 1+m2 = 
and the corresponding difference scheme u^j — ttj_ij_ mi — «j_ij_ m2 + «j_ 2 j- mi - m2 = 0. 
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For wave that containing N solitons from ^-presentations of the dispersion equation is 
possible to form following difference scheme. 



N M 



Yl a °n,i u i-nd-i = 0; (6) 



n=0 (=0 

a o,o = i; a o,i...M - °; a n, P - (-i) n ^ P ,5; 

n 

S = ^m CTi ; cr C m 1 . . .m N ; 
1=1 

n = l,2,...,N;p = l,2,...,M, 

where M is the sum of N velocity indexes m ar If p in (JHJ) is equal to the sum of n 
velocity indexes from given ensemble N indexes, then 5 - function is an unit, otherwise 
it's a zero. 

The double sum in expression (JBJ) presents itself equation of the two-dimention linear 
prediction (LP) jl] - samples of possible to define as linear combination of 

previous, given as matrix of size (iV+ 1) x (M+ 1). For iV solitons with arbitrary velocity 
indexes the LP model is given by 

N M 

^ ^ "„./''; i ~ 0; a 0i0 = 1, (7) 

n=0 1=0 

where M is rounded value of the sum of N velocity indexes. The LP parameters a n> i may 
be defined using the least square method on data of the field in initial layer. 

The derivatives on x and t can be presented as finite differences with use the equation 

©• 

1 / N 1 f , u it j , 51,0 



where 

AT M 
n=0 1=0 ( n ,l)^{0,0),(p,q) 

indexes in brackets denote the pair values. Such approach allows to increase the supporting 
area for calculation of the finite differences, as well as allows to normalize supporting areas 
of differences first and following orders. Using similar (JHJ), © presentations of the high 
order derivatives for the equation $1} shall receive following iteration difference scheme. 

Ui,j(k) = - (C oUi /(k - 1) + C 2 ) fd, (10) 
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where k=0,l,. 



- iteration step number; 



At f 1 



Oi j0 Ax z V a o,i ao,2 ai.o 01,2/ Ax \ a ,i 

a / 25*0,1 5*0,2 5"io , 25*11 , 5io , At / 5*0,1 0^0,1 

C 2 = "t 2 + + + + 7— 7 /^-j 2 - 

Ax V a ,i ao,2 «i,o ai,i 01,2 / «i,o Ax \ a ,i a^! 

In equation ([TU|) the nonlinear component is presented as discrete analogue of (3(u 2 ) x . As 
can be seen from ([10 J), the LP coefficients must be nonzero. For this the velocity indexes 
of the waves follow to choose as rational numbers with fractional part not multiple 2~ n . 

The numerical dynamic model of the initial layer allows to consider the interaction 
of solitons, moving in forward direction. The study of the models (J7J), (fTUj) has shown 
that the model (JJJ) reproduces the one-soliton wave with error linearly increasing from 
layer to layer, which relative value on a border of the net of size 500 x 500 not more 
0.001. The LP model reproduces two and more solitons with relative inaccuracy not more 
0.03, moreover, the inaccuracy value sharply increases in cross point of the solitons. Since 
given inaccuracy is a smooth function, reveals itself on solitons fronts and does't influence 
upon their form, that it possible consider as effect of the phase shift. The model (fT0|) 
reproduces the initial layer of the one-soliton wave with small phase shift. The iterative 
procedure was checked by convergence condition | Uij(k) — Uij(k — 1) |< l.e — 10 and 
was limited by ten steps. Performing the convergence condition was fixed in all nodes 
of the net. The studies have shown that important is the choice of the model order in 
accordance with (JHJ) - 2-image of (JZJ) must corresponds to the product of N dispersion 
equations (JHJ) each of solitons. Stability given model depends on magnitude of the LP 
coefficients, that depend on the choice of the sampling interval and velocity scaling. Is 
possible also the correction of the LP model order for reason to remove critical magnitude 
of the coefficients. 

On Fig.l is presented the interaction of two solitons of the RLW equation with 
parameters (0.5;0.125), when influence of nonlinearity is least. The velocities ratio less 
1.2. The model parameters are following: vq = 10; Aa; = 0.25; N = 3; M = 2; 
mi = 0.605; rri2 = 0.67. Before cross point solitons save the form, their amplitudes 
summarize linearly. Further occurs qualitative change of the wave form - are formed 
quick waves with slowly rising amplitude and a shock wave. Simultaneously with growing 
of the amplitude and narrowing the profile of the shock wave is formed a negative front, 
an amplitude which sharply increases after the positive front achieved a fixed level. The 
negative front moves with greater velocity and with increase the distance from saving form 
positive front it gradually forms a surge in area of great numbers. On correlation of the 
amplitude and velocity of the positive front, development of the negative front dynamic, 
smoothly moving over to singularity, possible draw the conclusion that the shock wave is 
transformed in a wave like (@J) accurate to phase shift and scaling multiplier on velocity 
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and amplitude. The numerical decision converges in accordance with presented above 
criterion in all nodes of the net. The integral value under profile of the wave was saved. 
Only in the field of sharply growing of the negative front the iterative process falls into 
closed cycle around neighbour values. If consider (fT0|) as quadratic equation then in this 
area its discriminant negative. Therefore the model in this area has evaluating nature. 
The form of singular area greatly is't changed with increase of quantity of iteration steps. 

Within the range of the velocities ratio 1.2. . . 1.25 the solitons interaction carries the 
elastic nature - the solitons disperse with conservation of the form and velocity. Given 
interaction possible to characterize and as a full mutual penetration, since in the cross 
point amplitudes of solitons summarize linearly. Similarly interact solitons with velocity 
indexes m2 = (K ± 0.25)mi, where K is integer value. 

On Fig. 2 is presented an instance of the solitons interaction with the velocities ratio 
outside of noted above areas: m\ = 0.55; = 0.7. The interaction occurs inelastically. On 
Fig. 3 is presented the consequence of shears of the wave before and after the interaction. 
From figures is seen that approximately only about half on amplitude and length of the 
quick soliton continues the motion with the previous velocity as the solitary wave. The 
length of this wave is a half of the length of the slow soliton. The slow soliton absorbed 
part of the quick soliton energy and transformed in the shock wave with the same scenario 
of the behaviour, as on Fig.l. 

If consider solitons as a sum of harmonical waves then from presented results follows 
that exist tunnel zones in neighborhood of multiple magnitudes of the velocities, when 
harmonic components mutually penetrate. Outside of the tunnel zones more slow wave 
emerges as potential barrier, through which passes the part of harmonical waves, period 
which is multiple with respect to length of the barrier. Rest harmonical waves enter in 
nonlinear interaction with the barrier and create the shock wave. Part of waves, which 
passed the barrier, creates stable on the form solitary wave, moving with previous velocity. 
As can be seen from the expression (J2J), such wave can be considered as solution of the 
equation (JTJ) with proportionally increased coefficients (a; (3). Consequently, resonance to 
the nonlinear medium can be SW, which parameters give appropriate balance of the form 
and velocity. 

In the equation with parameters (1;0.5), (l;-3) the nonlinear component has 
more deep influence. Therefore mutual penetration of the waves with small distortion 
was observed only in the tunnel zones with K > 2. On Fig. 4 is presented the solitons 
interaction of the RLW equation (1;0.5) with parameters of the model: N = 2; M = 3; 
rrii — 0.6; m-i = 0.7. The interaction runs with formation before front of the slow soliton 
the third wave, which gradually accumulates the main part of total energy of the solitons 
and develops in the shock wave, and next in wave like (@J). In this case the transition 
in the wave like (J1J) is realized quicker and so the shock wave excites the series of the 
fluctuations developing in periodical singularites. The third wave is formed as reflection 
of the slow wave from running up front of the quick wave. Then reflected wave absorbs 
quick soliton. Such interaction shows that with increase the influence of a nonlinearity 
grows the ability of the solitons to save the form. The interaction of the solitons of the 
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RLW equation with parameters (l;-3) occurs similarly and differs the inverse sign of the 
wave. 



Presented in article the model has allowed to see result of influence of the nonlinearity 
on non-equilibrium interaction of the SW. As a groundwork to influence is used the linear 
model (JZJ), which is equivalent to superposition of the operators (jHJ each of SW. 



where p n are integer values, which regulate the filtering characteristic each of multiplier. 
For realization of the computing scheme the auxiliary linear model can be formed 

for each layer or counting sample on its neighbourhood. The linear model must be unitary 
to not influence on processes of damping and excitation of fluctuations. With the help 
of linear model was eliminated ambiguity of nonlinear equation solution by matching 
on a dynamic of current and previous samples. Model was researched also with discrete 
analogue of the nonlinearity as — Wjj_2)/2A:c. Such model gave the results, 

similar mentioned above, but it was less stable in contrast with (fTUj) and amplitude of 
the shock wave did't rise up to saturation level. 

Possible suppose that transformation of the shock wave in the wave like (J1J occurs 
under influence of ju x in equation ((TJ. To value the level of this factor the model, 
presented on Fig.l, was researched for 7 = —0.25. . . 1.5. It was determined that ratio 
of the amplitude of the shock wave saturation to the average magnitude of the solitons 
amplitudes is valued as about 24—87 an d though amplitudes of the solitons with reduction 
7 increase the time of the shock wave saturation grows as about T{pf = 1) + 160A£(1 —7). 
These results show that 7 affects only partly - with the reduction of 7 decreases ability to 
saturation. Since a shock wave with high amplitude is't a solution of the RLW equation 
with any parameter 7, that the main factor is a process of self-organization of the shock 
wave under influence of the factors of amplitude, frontage and velocity, in consequence 
of which the shock wave shapes up the resonance form that is acceptable from condition 
of the high derivatives balance. This is confirmed following, when the solitons amplitude 
was changed by parameter t> = 10 . . . 200 the ratio of the amplitude of the shock wave 
saturation to average magnitude of the solitons amplitudes was valued as about 8t>o 0,3 
and the time of saturation greatly was't changed. However, was found ranges where does't 
occur the shock wave saturation: vq = 70 . . . 80; 125. . . 150; 180. . . 200. The shock waves 
with high amplitude are formed in these ranges right after collision, or in cross point of 
the SW. Whereupon the model lost stability. In addition, the double order was needed for 
stable linear model in these ranges. The model of the KdV equation behaved similarly. 

When v ^> 7 the spectral parameter in equation (jHJ) a ~ \/l/a and so solitons of the 
RLW equation are spatially quasicoherent. Hence, the nonequilibrium interaction of the 
quasicoherent solitons can generates solitary and shock waves, self-excited quick waves- 
harbingers. The shock wave can reach to be tenfold amplitude in contrast with generating 
it solitons due to the profile narrowing. In the case of weak nonlinearity possible partial 
mutual penetration with formation of saving form solitary waves. Strong nonlinearity 
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brings about forming the third wave, which accumulates the main part forming it the 
solitons. The self-excitation of the quick waves-harbingers can be considered as property 
of a nonlinear medium spontaneously to generate fluctuations under influence of small 
perturbations. 

The model has reflected conversion of the nonlinear medium from one type of the 
balance of the amplitude and velocity - two solitons (jHJ), to the other, in the form of 
the stable positive wave front, which precedes developing in singularity negative front. 
Physical sense of such transition possible to make clear considering other characteristic 
of concrete phenomena, described by equation ((Tj). For example, the RLW equation 
(0.5;0.125) describes an ion velocity in Broer - Sluijter's system of equations of plasma 
jS]. Numerical modeling of the dynamic of the ion density at the speed on Fig.l using the 
direct matrix method has shown that when ion velocity takes the form of the shock wave 
with stable amplitude and appears negative front the amplitude of fluctuations of the ion 
density sharply increases. The negative velocity possible to consider as effect of the elastic 
force reaction on quick ions. The rapprochement with the singular point of the velocity 
brings about forming the fluctuations of the density with high amplitude, spreading in 
forward and backward directions. Consequently, the singularity corresponds to excited by 
colliding ions source of the fluctuations of density. 

Solitary waves can meet in nature as multiwave systems and as separate formations, 
depending on this condition runs their interaction. The first type of the interaction well 
studied by analytical way. Example of the numerical analysis of the second type of the 
interaction is presented in this article. Results of modeling allow to draw a conclusion that 
nonlinearity reveals itself ambiguously - in one event exists the mutual penetration, in 
the other - full or partial interaction with transition in qualitative other resonance state. 
Partly the interaction of the SW in nonlinear medium possible to describe qualitatively as 
processes of the multiplying and mutual passing - reflections of their harmonic components 
taking into account orthogonality of the harmonics and conservations of total energy. 
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Figure 1: Solitons interaction with parameters of the model: a = 0.5; (3 = 0.125; 
Vx = 6.05; v 2 = 6.7; iV = 3; M = 2. 
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Figure 2: Solitons interaction with parameters of the model: a = 0.5; (3 = 0.125; v% = 5.5; 
r 2 7: A :>,: M = 2. 



8 



u(x,t) 

100 

5? SO 



O 




1 OO 200 



300 400 
grid t 



100 200 300 400 500 600 700 




Figure 3: Spatial shears of the solitons interaction with parameters of the model: a 
(3 = 0.125; vt = 5.5; v 2 = 7; N = 3; M — 2. 
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Figure 4: Solitons interaction with parameters of the model: a = 1; (3 = 0.5; v% 
r 2 7: N = 2; M = 3. 
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